Testing the Structure of Multipartite Entanglement with Bell inequalities 
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We show that the rich structure of multipartite entanglement can be tested following a device-independent 
approach. Specifically we present Bell inequalities for distinguishing between different types of multipartite 
entanglement, without placing any assumptions on the measurement devices used in the protocol, in contrast 
with usual entanglement witnesses. We first address the case of three qubits and present Bell inequalities that can 
be violated by W states but not by GHZ states, and vice versa. Next, we devise 'sub-correlation Bell inequalities' 
for any number of parties, which can provably not be violated by a broad class of multipartite entangled states 
(generalizations of GHZ states), but for which violations can be obtained for W states. Our results give insight 
into the nonlocality of W states. The simplicity and robustness of our tests make them appealing for experiments. 
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Entanglement |l|] plays a central role in quantum theory, 
in particular in the context of quantum information. Conse- 
quently, its characterization and detection are important is- 
sues, which have been extensively studied recently 0]. 

While the structure and detection of entanglement is rela- 
tively well understood in the bipartite case |l|,|2|], much less 
is known in the multipartite case (i.e. for three or more sys- 
tems). Given the importance of multipartite entangled states 
in the context of quantum computation yfl, quantum simu- 
lations quantum metrology @], and many -body systems 
0], an intense research effort has been recently witnessed in 
this direction. Sophisticated techniques, such as entanglement 
witnesses, have been developed (see e.g. J3l and references 
therein). This line of research is also important from an exper- 
imental perspective, given recent progress achieved with pho- 
tons fif], trapped ions jfj, and superconducting qubits ifioll . 

The central aspect of multipartite entanglement — which 
makes its study a highly nontrivial problem — is its rich and 
complex structure. Three (or more) quantum systems can be 
entangled in very different ways 111 ill . The most celebrated ex- 
ample is the case of three qubits, which feature two inequiva- 
lent classes of tripartite entangled states, namely Greenberger- 
Horne-Zeilinger (GHZ) Q2] and W states Ultl . 

Importantly, different types of multipartite entanglement 
can have completely different features, for instance from the 
point of view of strength of correlations, or robustness to 
noise and losses. In general this greatly affects their informa- 
tion theoretic properties. For instance, GHZ states, which are 
highly entangled but very fragile with respect to losses, are es- 
sential for many quantum information-theoretic tasks, such as 
quantum secret sharing lfl3ll and quantum teleportation 1 14]. 
On the other hand, W states, which are much less entangled 
but highly robust against losses, are central to the physics of 
the interaction of light and matter 1151 . in particular for quan- 
tum memories 11611 . More generally, it is also known that too 
much entanglement hinders computational power Il7ll . 

It is therefore important to understand how different types 
of multipartite entanglement can be characterized and distin- 
guished from each other. Recently Schmid et al. Ill 811 and 
Huber et al. Il9ll discussed this problem, introducing entan- 
glement witnesses that allow one to distinguish between im- 



portant classes of qubit states. 

However, entanglement witnesses have an important draw- 
back, in particular from the experimental point of view, which 
is that their derivation is based on assumptions which are dif- 
ficult (if not impossible) to meet in practice. In particular, the 
violation of an entanglement witness can be considered as a 
conclusive proof of entanglement only if the measurement op- 
erators are precisely characterized. This is indeed impossible 
in practice where (small) alignment errors are inevitable. Im- 
portantly, such errors can have dramatic consequences, as re- 
cently shown in Refs I20I1 ; uncertainties on the measurements 
make it possible (in principle) for a separable state to lead to 
a violation of the witness, thus rendering the test inconclu- 
sive. Therefore, it is desirable to devise methods for both wit- 
nessing multipartite entanglement and for testing its structure, 
which do not rely on any assumptions about the measuring 
devices used in the protocol. 

The first problem has been addressed recently via the intro- 
duction of device-independent witnesses for multipartite en- 
tanglement Ir20ll . Here, we address the second question, and 
show how the structure of multipartite entanglement can be 
tested in a device-independent manner. To do so, we present 
Bell inequalities which (i) are tailored for detecting the non- 
locality of certain multipartite entangled states, and (ii) can 
provably not be violated by other classes of multipartite en- 
tangled states. Hence the violation of our inequalities allows 
one to discriminate between various types of multipartite en- 
tanglement. Importantly this is achieved without placing any 
assumptions on the measurement devices used in the protocol, 
since the violation of a Bell inequality is a statement based 
only on the observed statistics. Here, we start by address- 
ing the case of three qubits, and show how one can distin- 
guish between GHZ and W entanglement. Then we move to 
the more general case and introduce a class of Bell inequal- 
ities which can provably not be violated by a large class of 
multipartite entangled states (in particular generalized GHZ 
states). The key feature of our inequalities is that they are 
based on sub-correlations, i.e. each term of the Bell polyno- 
mial involves only a subset of the parties. We derive a family 
of 'sub-correlation Bell inequalities' for an arbitrary number 
of parties, and show that they can be violated by W states. Fi- 
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nally we describe a method for deriving useful sub-correlation 
Bell inequalities and apply it to simple cases. Generally, our 
results give insight into the nonlocality of W states [21], which 
is still poorly understood. 

Preliminaries. We will focus on two inequivalent classes 
of multipartite entangled states of N systems: (i) W states: 



1 



,01) + |0... 10) + . .. + 110...0)) (1) 
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Here, following lfl9ll . entanglement classes are defined by 
equivalence under local unitaries and permutation of the par- 
ties 13 ill - Both classes are robust to imperfections, however in 
a very different manner. For instance, GHZ states are maxi- 
mally entangled and thus robust against noise. On the other 
hand, W states are robust against losses; tracing out k <C N 
particles leaves the state barely unchanged. This is not the 
case for generalized GHZ states which are extremely fragile 
to losses; indeed it suffices to trace out a single subsystem to 
destroy all entanglement. This difference will play a crucial 
role later on. 

The central tool here will be Bell-type inequalities. We 
shall focus on multipartite Bell polynomials which are sym- 
metric under any permutation of the parties. Each party can 
perform two possible measurements featuring binary outputs. 
Therefore the following notation will be useful; we write 
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a Q (Ao + B ) +a 1 (A 1 + B x ) 
+a aQ A a B Q + a i(A Si + AiB ) 
+a 1 iA 1 B 1 (4) 

where Ai = ±1 denotes the outcome of measurement setting 
i = 0, 1 of party A. Likewise for Bi. The extension to more 
parties is straightforward. For instance, for N = 3 parties, the 
Mermin inequality M22I1 . usually given by 

M 3 = A B C - A B 1 C 1 - A x BqC x - A x B x C a < 2 (5) 
now reads 

M 3 = [0 ; ; 1 - 1 0] < 2. (6) 
Similarly Svetlichny's inequality l23ll will now be written as 

S 3 = [0 0;000;ll-1-1]<4. (7) 

while the generalized Mermin-Ardehali-Belinskii-Klyshko 
[Hill] (MABK) Bell inequality is given by 
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Case of 3 qubits. We show how to discriminate in a device- 
independent way between the states \ W%) and IGHZ3). 

First, we consider the Mermin inequality ©. Here the 
state IGHZ3) can lead to maximal violation (i.e. A/3 = 4), 
which corresponds to the GHZ paradox lfl2Tl . However, only 
much lower violations are possible for the |W 3 ) state. Fol- 
lowing the method of l26ll . based on semide finite program- 
ming, we could prove that A/3 < 3.079 for | W3); details can 
be found in Appendix A. Thus a violation of the inequality 
A/3 < 3.079, which is possible using IGHZ3), certifies that 
no \Ws) states were used. The robustness of this effect is 
best represented using the resistance to noise, which is here 
found to be wghz 3 ~ 77%. Physically this means that unless 
a fraction of (at least) (1 — u>ghz 3 ) of white noise is mixed to 
|GHZ 3 ) the inequality A/ 3 < 3.079 will be violated. 

Next, we study the converse task. Let us consider the fol- 
lowing inequality 

B = S 3 + R 2 = [ ; 1 ; 1 1 - 1 - 1] < 6. (9) 

Thus B is simply the Svetlichny polynomial S3 ©, to which 
we have added some 2-party correlation terms, denoted i?2- 

On the one hand, it turns out that inequality (|9]l can be vi- 
olated by the state IW3). The largest violation is B ~ 7.26 
(see Appendix B), leading to a resistance to noise of Ww 3 — 
82.7%. On the other hand, inequality (O can provably not be 
violated by generalized GHZ states (0 with d = 2; we con- 
jecture that inequality (0 holds for all states (0, which we 
checked numerically for d < 10. Below we sketch the proof 
for the state IGHZ3); the full proof for 2-dimensional gener- 
alized GHZ states can be found in Appendix C. Finally, note 
that inequality (|9]i holds for any bipartition, that is, for any 
configuration in which two (out of three) parties can group 
and have any no-signaling correlations. 

Sketch of the Proof. Consider the decomposition of the 
polynomial B = S3 + /?2- The intuition behind the proof 
is the following. On the one hand, the state IGHZ3) can lead 
to a hig h violation of the Svetlichny inequality (S3 < 4\/2, 
see Il25l0 using qubit measurements in the XY plane (i.e. the 
equator) of the Bloch sphere. In that case however, the state 
will exhibit no 2-party correlations, leading to R2 = 0. Thus 
overall we get no violation. On the other hand, the GHZ state 
has perfect 2-party correlations for measurements along the Z 
axis. Thus one can achieve R2 = 6 with such measurements. 
However, in this case, there will be no 3-party correlations, 
thus leading to S3 = 0. Overall, we get no violation again. 

Based on this intuition we build the proof as follows. First 
we write the measurement operators as Aj = cos 9 jXY + 
sin 6j Z where XY denotes a Pauli operator on the XY plane, 
i.e. of the form cos 4>X + sin <j)Y\ here X, Y, Z are the usual 
Pauli matrices. One can then show that the S3 part of the poly- 
nomial is upper bounded by an expression depending upon the 
cosine of the angles 6j, whereas the R2 part of the polynomial 
is upper bounded by an expression of the sine of the angles 9j . 
After trigonometric massaging, it follows that B < 6. ■ 

General case; Sub-correlation Bell inequalities. We now 
move to the case of an arbitrary number of parties. Our main 
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tool will be sub-correlations Bell inequalities (SCBI). These 
are based on polynomials which do not include any full N- 
party correlation terms. Here we focus on SCBI which are 
symmetric (under permutation of the parties). 

The motivation for studying such inequalities comes from 
the fact that different types of multipartite entanglement are 
affected differently by losses of particles. On the one hand, 
losses barely affect W states. For N large enough, we have 
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where pw N = \ Wn)(Wn\, which suggests that W states can 
violate SCBI. On the other hand, one would intuitively ex- 
pect that states which are not entangled anymore after having 
traced out one particle, will never violate a sub-correlation 
Bell inequality. Although this may not be true in general 12711 . 
we will see below that it holds nevertheless for a wide class of 
multipartite entangled states, in particular for all generalized 
GHZ states ©. 

In order to show that a given state p^ does not violate any 
SCBI, it is sufficient to ensure that all of its N reduced states, 
obtained by tracing out one party, are compatible with a sepa- 
rable state <7jv, i.e. \ik{PN) = trfc(djv) for all k G {1, N}. 
In this case, there can be no Bell violation, since the separable 
state it n provides a local model. For any state of the form 
@, a separable state satisfying the above condition is given 

by <Jn = Ej=o I a 3 \ 2 1 3 ■ ■ ■ j) ti ■ ■ ■ 3 1 ■ Thu s generalized GHZ 
states can never violate any SCBI. 

We now present a family of sub-correlation Bell inequal- 
ities for an arbitrary number of parties N, which can be vi- 
olated by W states (UJ). This construction is based on the 
MABK polynomials, and is given by 
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where M^_j denotes the polynomial in ([8]l involving all par- 
ties except party k. Given that M^_ a < 2^^^! (see ©), it 

follows that Bn is upper bounded by N2^^~\ the sum of 
all the local bounds. It turns out that this bound is tight, i.e. 
there exists a local model (for all N parties) which achieves 
the local bound of each MABK polynomial in (fTTT >. 

Next we study the achievable violations of inequality (fTTT i 
for states |Wjv}' We start by writing the W density matrix in 
the basis of Pauli matrices { 1 , X, Y, Z}\ 



TABLE I: Violations (given as resistance to white noise) for | Wn) 
states of the SCBI flTTl . 
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10 12 15 20 40 



1 0.891 0.831 0.792 0.765 0.730 0.709 0.688 0.669 0.642 
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+2(l + Z) CSN - 2 (XX + YY) 



where sym[a;] means the symmetrization of expression x. 

Noting the XY symmetry in the above expression, we focus 
here on measurements in the XZ plane of the Bloch sphere 
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0, 1). Also, we assume that all parties perform the same two 
measurements, i.e. Aj = Bj = Cj and so on. 

We can now derive an expression for a full correlation term 

A x (8> By <g> . . . ® N z featuring N — k settings '0' and 
k settings ' 1 ' : 
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Inequalities ( fTTT i involve only correlation terms of (JV — 1) 
parties. Using the equality in ( [Tol l, we can evaluate a term 
1 — k settings '0' and k settings T: 



T;^ 1 featuring N 



1 

N 
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(14) 
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The first line of the above expression is then given by 
Cq dy/N, while the second line can be evaluated using ( TPl ). 

Thus, we get an expression for the value of B^ for states 
|Wjv) which is a function of only two angles (^o,i)- Using 
numerical optimization, we investigated W violations up to 
TV = 40. The results are presented in Table 1, and are pre- 
sented in terms of resistance to noise, which is given here by 
ww N = A^2 1" /tr(pw N BN). The violation of dTTb thus 
clearly increases with N; we conjecture that for all > 5, 
the W state leads to a Bell violation. Also, for N = 5, we 
could check numerically that our choice of measurement set- 
tings is optimal. 

Finding sub-correlation Bell inequalities. Further relevant 
SCBI can be found using the geometrical approach to non- 
local correlations. In the latter, Bell inequalities correspond 
to the facets of the polytope of local correlations. The ex- 
tremal points of this polytope are given by local deterministic 
strategies — in which each party outputs deterministically for 
all measurements — and are thus easily characterized. Then, 
using appropriate software, it is possible to find the facets of a 
given polytope, and thus retrieve Bell inequalities. 

In order to find SCBI, we follow the above approach, but 
forget about full (A-party) correlation terms, thus effectively 
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projecting the initial local polytope in a subspace of smaller 
dimension. Here we have followed this method in the sym- 
metric approach (i.e. focusing on inequalities which are sym- 
metric under permutation of the parties), which greatly sim- 
plifies the problem lf28ll . and considering two measurements 
per party, all with binary outcomes. This allowed us to find 
SCBI for N = 4, 5. Here we present two of them, which are 
relevant from the point of view of quantum violations: 



h = [-1 -l;-2 -2; -2 11 -2]<8 (15) 
h = [0 0; -2 -1;0000;-4 0201]<15 

Note that here we have used the notation (@}, omitting the co- 
efficients for full correlation terms since the latter are always 
zero for SCBI. For these inequalities, we have found large vi- 
olations for the W state; for I4 we get wwi — 0.707, while 
for I5 we get ww s — 0.536 (see Appendix B). 

Finally, it is worth mentioning that a phenomenon of frus- 
tration occurs for both of these inequalities. Specifically, each 
symmetric SCBI Jjy < L can be viewed as the symmetriza- 
tion (over N parties) of a Bell inequality jV-i < I involving 
only N — 1 parties (see eq. (fTTTi). It is then natural to define 
a frustration parameter F = ^P. When F = 1, there is no 
frustration, since there is a local deterministic model which 
saturates the local bound of each inequality (i.e. Jn-x = 
which thus reaches the overall upper bound of Jjy = L = Nl. 
This occurs for instance in our construction of eq. ill) . In- 
terestingly, this is not the case for the two inequalities above 
(fTBT l, which feature a rather large frustration: F = 11/3 in 
both cases. It would be interesting to see whether large frus- 
trations are favorable to large quantum violations. 

Conclusion and Perspectives. We have addressed the ques- 
tion of distinguishing between different classes of multipar- 
tite quantum entanglement in a device-independent manner. 
In particular, we have presented Bell inequalities which can 
provably not be violated by generalized GHZ states but can 
be violated by W states for an arbitrary number of parties. We 
believe that the simplicity and robustness of these tests make 
them attractive from an experimental viewpoint. 

Finally let us comment on some perspectives for further 
work. First, for N > 4, there exist many other classes of en- 
tangled states I29I1 . such as Dicke states fl5ll and cluster states 
J^t], for which it would be interesting to find similar tests. 

Also, we believe that these ideas, in particular the concept 
of sub-correlation Bell inequality, may lead to further insight 
into multipartite quantum nonlocality. For instance, they may 
allow one to better understand the nonlocal correlations of the 
cluster state l30ll . which could give insight to the computa- 
tional power of the latter which is still not fully understood. 
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APPENDIX A 

Here we derive an upper bound on the violation of the Mer- 
min inequality M3 < 2 with the state | W3). Using the method 
described in Ref. 12611 . based on semidefinite programming 
[see e.g. J. Lasserre, SIAM J. Optim. 11, 796 (2001)], we 
prove that M 3 < 3.0792 for \W 3 ). 

We start by writing the measurement operators in the form: 



= ajX + faY + jjZ 



(16) 



where on 



ft 



1. Note here, since we consider a 
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Bell inequality featuring binary outcomes and 3-qubit state, it 
is enough to consider qubit projective measurements. 

Next, we write the density matrix p\y 3 = | W3) 1 in the 
basis of Pauli matrices, i.e. {t,X,Y,Z}. We find: 



PWs = 24(3111 + Zti - ZZt 

- 3ZZZ + 2(1 + Z)(XX + YY)) 



(17) 



where we have omitted symmetric terms. 

Our task now is compute an upper bound on the expression 
E(u) = tT(Mspw 3 ), where M.3 denotes the Mermin Bell 
operator, i.e. the operator we obtain from inserting measure- 
ment operators of the form ( TTBT ) into the Mermin polynomial 
(see main document, eq. (5)), and u denotes all the parameters 
characterizing the measurement operators. Thus, our problem 
is to solve 



E* = max u E(u) s.t. c(u) > 



(18) 



where E(u) is qubic function of u, featuring quadratic 
constraints on it, denoted c{u). This problem can 
be solved using semidefinite programing. Here we 
have used the package GloptiPoly3 [available at: 
http ://homepages . laas . fr/henrion/software/gloptipoly 3/] . 

We obtained E* < 3.0792 in the second order of relax- 
ation. This bound is most likely not tight. Numerically 
we have found that the optimal violation for the state \Ws) 
is M3 ~ 3.0460, which leads to a resistance to noise of 
ww 3 — 0.6566. The measurements can be chosen to be real 
(i.e. in the XZ plane of the Bloch sphere) and the same for all 
parties, that is, of the form 



A; = B, 



cos6jZ 



iX. 



(19) 



The optimal angles are given by 9q = 0.30027T and 9\ = 
0.8673tt. 

As it was mentioned in the main text, the Mermin in- 
equality is maximally violated by the GHZ state. The op- 
timal measurements are those featured in the GHZ, namely 

A = Bo = Co = X and A 1 = B l = d = Y. 



Here we give more details about the quantum violations, 
in particular for W states, of all the relevant inequalities pre- 
sented here. In particular, we give the optimal measurement 
settings, which can always be taken to be real (i.e. in the XZ 
plane of the Bloch sphere), of the form 



Aj = cosQjZ 



sniOjX. 



(20) 



Moreover, it is also optimal to consider identical measure- 
ments for all parties, i.e. Aj = Bj = Cj and so on. 

First, we consider the Bell inequality B = S3 + R2 < 6. 
The optimal violation for the state IW3) is B ~ 7.2593, 
which leads to a resistance to noise of ww 3 — 0.8265; 
the measurements are given by 6q = 0.2677-7T and 9% = 
n — 60 ■ Note that the maximal violation for qubit states is 
slightly larger. We have found B ~ 7.3084 for the state 
|V>) = 0.99711^3) - 0.07597|111) by measurement angles 
9o = 0.2615tt and 6>i = n - 6 . 

Next we move to the sub-correlation Bell inequalities (see 
main document, eq. (15)). For inequality I4 < 8 the largest 
violations for the state \W&) is found to be B ~ 11.3155, 
which leads to a resistance to noise of ww 4 — 0.7070; the 
measurement settings are given by 60 = 0.78617T and 9\ = 
2tt — 9o- The true quantum maximum is B ~ 12.0680, which 
occurs for the state \tp) = 0.9877|W 4 ) - 0.1561|W4) by the 
following measurement angles #0 = 0.76657T and 6>i = 27r — 
6*o, where \W4) is the spin-flipped | W4) state. 

For inequality 7g < 15, the largest violations for the state 
IW5) is found to be B = 28, which leads to a resistance to 
noise of ww 5 — 0.5357; the measurement settings are given 
by 6>o = and 9\ = it/ 2. On the other hand, the true quan- 
tum maximum is B ~ 30.1918, which occurs for the state 
\ip) = 0.97528|Ws) - 0.2201|D|) - 0.01851|11111) by the 
measurement angles, 9$ — and 9\ = tt/2, where |£>|) is the 
5 -qubit Dicke state with three excitations. 



APPENDIX C 

In this Appendix, we give the proof that the Bell inequality 
B = S3 + i?2 < 6 (see main document, eq. (9)) cannot be 
violated by generalized GHZ states (see main document, eq. 
(2)) of dimension 2. 

Proof. We start by proving the above bound for the 3-qubit 
GHZ, which can then be extended to all 2-dimensional gener- 
alized GHZ states at the end of the calculation. 

We first write the density matrix for the state in terms of the 
Pauli matrices. Using the relations 



|0)(0| = 
|0)(1| 



l + Z 
2 

X + iY 



|1)(1| = 
U)(0| - 



1 - z 

2 

X-iY 



(21) 



(22) 
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we have that 



p=-(|000) + |lll))((000| + (lll|) 

= -{ttt + (ZZt + ZtZ +1ZZ) (23) 
8 

+ XXX - (XYY + YXY + YYX)) 

Next, we wish to parameterise our measurement operators 
that appear in the polynomial B. The mth measurement on 
the qubit j will be parameterised as: 



M] n = cfXY + sfZ 



(24) 



where cj 1 = cos(i™), s™ = sin(i™) and XY denotes a Pauli 
operator on the XY plane. I.e. of the form cos <fiX + sin <f>Y. 

Using this parameterisation and substituting in to tr(pB) 
yields the value 



ti(pB) = c?c°c° + c?c°c* + c?4 C o + c\clcl 

— 444 — 4c 2 4 — C 1 C 2 C 3 — 444 
+ s l4 + 4 S 2 + s l s 3 + 44 + 44 + 44 



(25) 



We now wish to show that this expression is upper bounded 
by 6. We will first optimise over t\ and t\, which may be done 
explicitly using the identity 



a sin 9 + b cos 8 — \/ a 2 + b 2 sin(# + i 



where 



arctan ( - ) 



if a > 



[ arctan (^) + n, if a < 
so the maximal value of this must be 



(26) 



(27) 



maxe{a sin 9 + b cos 9} = V a 2 + b 2 (28) 
Hence we can write 



tz{pB) < y/a 2 + b 2 + V ' c 2 + d 2 + s§4 + 5^° (29) 
where 



a = S n + Si 



b — c 2 c 3 + c 2 c 3 + c 2 c 3 c 2 c 3 



^ — C 2 C 3 C 2 C 3 C 2 C 3 C 2 C 3 



(30) 
(31) 
(32) 
(33) 



As s 2 s 3 + s 2 Sg has a maximum algebraic value of 2, it suffices 
to show that 



Va 2 + b 2 + \Jc 2 + d 2 < 4 (34) 
or equivalently (by squaring the previous equation) 



a 2 + b 2 + c 2 + d 2 + 2 v / (a 2 +b 2 )(c 2 + d 2 ) < 16 (35) 



We will first show that 

a 2 + b 2 + c 2 + d 2 < 8 (36) 
which can be expanded to 

a 2 + b 2 + c 2 + d 2 < 4 + (c° 2 c° 3 ) 2 + {c\c\) 2 
+ (c 2 , 4) 2 + + 244 + 244 < 8 (37) 

meaning that we only have to show 

(44? + (4clf + (44) 2 + (44) 2 

+2^4 + 2s^4 < 4 (38) 

In order to do this we will find all stationary points of the 
function and show that at these points the function does not 
exceed 4. Taking half the terms and maximising over £ 2 yields 

max t o{(c 2 4) 2 + (c°4) 2 + 244} 

= ((4) 2 + i)(4) 2 + 244 (39) 

Differentiating this statement with respect to t\ and setting it 
equal to zero allows us to find 2 conditions for a stationary 
point, either 



or 



4 = 



s 3 ~ s 2( c 3) 



(40) 



(41) 



If c 2 = then equation ( f39l > simplifies to 1 ± 2s\ — (s\) 2 , 
which takes a maximum value when s\ = ±1 and hence equa- 
tion ((39]) < 2. On the other hand if S3 = s 2 (c3) 2 then equation 
( |39l becomes 2(03) 2 + tan 2 (t|), which has a stationary point 
only when s| = and hence the maximum must be 2. 

The other half of equation (l38l takes the same form and is 
proved in the same way. Therefore 



(44) 2 + (44) 2 + (44? + {44? 

+244 + 2 44 ^ 4 

All that remains to be done is to show that 



2vV + b 2 )(c 2 + d 2 ) < 8 
or equivalenty (again by squaring) 

(a 2 +b 2 )(c 2 + d 2 ) < 16 



(42) 



(43) 



(44) 



In order to show this we will make use of the previously 
shown fact that a 2 + b 2 + c 2 + d 2 < 8. Let / = a 2 + b 2 and 
g = c 2 + d 2 , then we need to show that fg<16 given that 
f + g < 8. This last equation can be rearranged to find 



9<8-f^fg<f(8-f) 



(45) 



which is maximised when / = 4, which means that the max- 
imum value of fg must be 16. This final part proves that 
tr(pB) < 6 for GHZ states. 
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We now extend this result to all generalized GHZ states of 
dimension 2: 



p = M000) + flflllll)) (c e (000| + 

= - [111 + {ZZt + ZIZ+ IZZ) 
8 

+ c 2 e{tZl +Z11 + ttZ + ZZZ) 
+s 2e (XXX - (XYY + YXY + YYX))} 



(46) 



It is possible to maximise this statement over analytically 
by finding the function's stationary points. Differentiating and 
setting it equal to zero means either 



or 



C 1 C 2 C 3 



" C 1 C 2 C 3 



C 1 C 2 C 3 



C 1 C 2 C 3 



C 1 C 2 C 3 



C 1 C 2 C 3 



C 1 C 2 C 3 



c\c\c\ 



(48) 



Therefore we can calculate the value of tr(pB) using our pre- 
vious measurement parameterisation 



tv(pB) = s 2S (c { l4c u 3 + 



C 1 C 2 C 3 



C 1 C 2 C 3 



1„0 „0 



1 L 2°3 



C 1 C 2 C 3 



,1X1 
-1 L 2 C 3 



r l r l„0 
C 1 C 2 C 3 



s 1 s 2 + s 1 s 2 + s 1 s 3 



C 1 C 2 C 3) 



'2 a 3 



(47) 



'2 a 3 



cos(26») = 



(49) 



If equation (l48b is true then the equation reduces to a sum of 
6 values, each less than or equal to 1, so the bound holds. If 
equation (l49l holds then 9 = f and the state reduces to a GHZ 
state, where we have already proven the bound. Therefore this 
bound of 6 holds for all Schmidt states as well. 



